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Abstract 


Vibration and noise caused by gear dynamics at the meshing teeth propagate through 
power transmission components to the surrounding environment. The purpose of this 
work is to develop computational tools to investigate the vibro-acoustic propagation of 
gear vibration and to investigate the effect of different bearing types on noise radia- 
tion. Detailed finite element /contact mechanics and boundary element models of the 
gear/bearing/housing system are established to compute the system vibration and noise 
propagation. Both vibration and acoustic models are validated by experiments including 
vibration modal testing and sound field measurements. 

Bearings are critical components in drivetrains. Accurate modeling of rolling element 
bearings is essential to assess vibration and noise of drivetrain systems. This study also 
seeks to fully describe the vibro-acoustic propagation of gear dynamics through a power- 
transmission system using rolling element and fluid film wave bearings. Fluid him wave 
bearings have higher damping than rolling element bearings and so could offer an energy 
dissipation mechanism that reduces the gearbox noise. The effectiveness of each bearing 
type to disrupt vibration propagation is explored using multi-body computational models. 
These models take into account gears, shafts, rolling element and fluid him wave bearings, 
and the housing. Radiated noise is mapped from the gearbox surface to surrounding 
environment. The effectiveness of each bearing type to disrupt vibration propagation is 
speed-dependent. Housing plays an important role in noise radiation. It, however, has 
limited effects on gear dynamics. 

1 Introduction 

Gearbox vibration contributes the structural-borne noise in helicopters [9]. The gearbox 
noise consists of a wide range of gear mesh, shaft, and bearing frequencies within relative 
lower audio frequency range than the jet engine noise, which is another source of helicopter 
structural-borne noise. Systematic studies of the noise and vibration behavior are essential 
to design quiet helicopters and to identify noise and vibration sources. Limited work has, 
however, investigated the relationship between the gearbox noise and vibration. 

Experimental vibro-acoustic analysis of helicopters requires effort [9, 10, 11, 12, 13] 
because advanced signal processing techniques are needed to separate the gear, shaft, and 
bearing signals [14]. The analytical or computational vibro-acoustic analysis of geared sys- 
tems is also sparse [15] in the literature due to the complexity of these problems. Structural- 
borne noise calculations of geared systems are, furthermore, semi-empirical [16, 17]. These 
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acoustic models are not able to capture realistic dimensions of gearboxes which often 
have complicated structures, which affect the noise estimate. This work defines the vibro- 
acoustic behavior of gear dynamics through power-transmission systems using multi-body 
dynamics gearbox models established in [18, 19] and an in-house program. 

Dynamic forces at the meshing teeth drive a system vibration through power trans- 
mission components to the fuselage. The bearings linking the gear shafts to the housing 
are a primary factor in this noise path. Fluid film wave bearings are a special type of 
journal bearings, which have waved inner diameters of the stationary bearing sides shown 
in Figure 1. The wave profile amplitude is a few micrometer so that it is not obvious 
in this figure. The wave bearing technology used for gas turbine lubrication [20, 21] is 
recently applied to the planet bearings used in aviation planetary gears [22, 23, 20]. This 
technology provides higher stiffness and better lubrication for the bearings. Experiments 
on an aviation gearbox [22] with wave bearings demonstrate 25% higher load capacity com- 
pared to plain journal bearings. Dimofte [20] compared the load capacity between wave 
and journal bearings through an analytical formulation. He concluded that wave bearings 
are more stable than plain journal bearings under light-load or unloaded conditions in any 
operating regime. Furthermore, wave amplitude and starting positions of the wave profiles 
are important parameters affecting wave bearing performance. 

Machinery applications use rolling element bearings that do not create meaningful 
damping to reduce the transmitted structural-borne noise. Wave bearings have higher 
damping [24] and could offer energy dissipation. If fluid Him wave bearings are shown to 
withstand harsh operating conditions and provide better vibration characteristics, they 
may prove an attractive alternative to standard rolling element bearings. 

The major objectives of this study are to: 1) develop a finite element gearbox model 
which includes the detailed contact analysis of the gear tooth mesh and individual bearing 
rolling elements; 2) build up analytical (lumped parameter) model of the gear/bearing/housing 
system, which provides efficient dynamic analysis; 3) establish a boundary element model 
of the gearbox housing and map the radiated noise from the gearbox surface through 
acoustic analysis; 4) validate the vibration and acoustic models of the examined gearbox 
against measurements and theoretical solutions in the literature; and 5) understand the 
effectiveness of fluid film wave bearings in breaking the vibro-acoustic propagation path 
from the gears to the housing. 
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Figure 1: Fluid film wave bearing. 

2 Preliminary Step: Experimental Modal Test 

Frequency response measurements provide vital system information: natural frequencies, 
mode shapes, and damping characteristics. Modal analysis by impact testing establishes 
a baseline of comparison for mathematical simulations. Computational and analytical 
tools are used to fully characterize the performance of rolling element and fluid film wave 
bearings. Experimental measurements provide benchmarks, against which mathematical 
models are validated. 

Impulse test data was collected on the NASA GRC gear test rig by Ohio State re- 
searchers in November, 2008. Data was collected to provide frequency response curves of 
the gear test rig for multiple inputs and outputs. This data provides transfer functions 
across system components and can be used to extract vibration mode shapes. A stan- 
dard impact hammer with built-in force transducer was used to provide a known input 
excitation. Accelerometers were used to obtain dynamic response. The known input and 
measured output yield system frequency response characteristics that are independent of 
the testing method. These results help predict the system behavior in response to gear 
vibration, establishing benchmarks that validate the modeling effort. The experimental 
study also helps defining the system for computational/analytical study. 

At the time of the OSU experiments, four accelerometers were already installed by 
NASA researchers. One accelerometer was mounted vertically and one horizontally on 
each of the pinion shaft bearings. Nine additional accelerometers were mounted and used 
by OSU to measure the frequency response around the test gearbox. Figure 2 shows the 
approximate locations of these accelerometers mounted in and around the gearbox. Impact 
tests were performed throughout the system. This report contains findings from four of 
the most relevant impact tests. Figure 3 identifies the impact locations for these four tests. 
Three impacts were averaged for each experiment. Data was taken with the gearbox cover 
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off to measure vibrations inside the gearbox. The subsequent results are shown on a dB 
scale. 



Figure 2: Nine mounting locations for additional OSU accelerometers in and around gearbox. 

2.1 Translational Impact to Gear Tooth 

In this test an impact was made to a gear tooth to simulate excitation near the meshing 
gears. Figure 3(a) shows the impact point and two relevant accelerometers mounted on 
the pinion shaft (A5 and A6). The averaged frequency response function from the NASA 
accelerometers is shown in Figure 4(a). The frequency response functions from the relevant 
OSU accelerometers for this test are shown in Figure 4(b). 

The data in Figure 4 shows distinct natural frequencies at 750 Hz and 2500 Hz in the 
NASA bearing accelerometers (a) and the additional accelerometers (A5 and A6) mounted 
to the pinion shaft (b). The pinion shaft also shows a smaller peak around 1800 Hz. Both 
graphs show less-defined dynamic behavior between 3300 Hz and 3800 Hz. Accelerometers 
A1 and A2, which are mounted to the input shaft pillow block outside the gearbox (Figure 
2), show a reduction of vibration amplitudes by about 20 dB. This shows that vibration 
propagation is well-contained within the gearbox. Most importantly, this test shows two 
natural frequencies associated with the gears and/or shafts at 750 Hz and 2500 Hz which 
appear to be contained within the gearbox. 
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(a) 


(b) 




(c) (d) 

Figure 3: Impact locations for relevant tests: (a) translational impact to gear tooth with ac- 
celerometers A5 and A6 mounted to the pinion shaft, (b) impact to output gear shaft with 
accelerometers A5 and A6 mounted to the pinion shaft, (c) impact to rig input shaft with ac- 
celerometers A8 and A9 mounted to the input shaft, (d) impact to rear of gearbox housing with 
accelerometers A3, A4 and A7 mounted to the gearbox housing. 

2.2 Translational Impact to Gear Shaft 

In this test an impact was made to the output gear shaft, between the gear and the 
supporting bearings. This impact location is near the previous test, and it was performed 
to determine if the dominant 750 Hz and 2500 Hz peaks are predominantly gear or shaft 
modes. Figure 3(b) shows the impact point and two relevant accelerometers mounted on 
the pinion shaft (A5 and A6). The averaged frequency response function from the NASA 
accelerometers and the frequency response functions from the relevant OSU accelerometers 
for this test are shown in Figure 5(b) and compared to the results from the previous test 
to the gear tooth in Figure 5(a). 

The results from the test on the gear shaft show the same two dominant peaks at 
750 Hz and 2500 Hz. The amplitudes, however, have changed significantly. The 750 Hz 
peak is about 5 dB higher in the impact to the gear tooth, but the 2500 Hz peak is about 
10 dB higher in the impact to the gear shaft. This suggests that the 2500 Hz mode is 
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Figure 4: Frequency response functions for impact test to output gear tooth with accelerometers 
A5 and A6 mounted to the pinion shaft and accelerometers Al and A2 mounted to the input 
shaft pillow block for (a) mean NASA accelerometers and (b) select OSU accelerometers. 


dominated by shaft vibration and the 750 Hz mode has more gear body translation. 
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s' -10- 
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Frequency (Hz) Frequency (Hz) 

(a) (b) 

Figure 5: Frequency response comparison for (a) impact to pinion gear tooth, and (b) impact 
to gear shaft. 





2.3 Impact to Gearbox Input Shaft 

In this test an impact was made to the large, hollow input shaft of the test rig. This test was 
performed to measure vibration transmissibility into and out of the gearbox. Knowing that 
this shaft (and the output shaft like it) will have bending and torsional modes within the 
frequency range of interest, it was necessary to determine if these accessory components 
were coupled to the gearbox system under study. Figure 3(c) shows the impact point 
and two relevant accelerometers mounted on the input shaft itself (A8 and A9). All four 
frequency response functions from the NASA accelerometers are shown in Figure 6(a). The 
frequency response functions from the relevant OSU accelerometers for this test are shown 
in Figure 6(b). 
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The data in Figure 6(b) shows five distinct natural frequencies of the input shaft be- 
tween 3050 H z and 4500 Hz, measured by A8 and A9. This is within the frequency range of 
expected gear dynamics. Therefore, if the input shaft is coupled with the gearbox system, 
it would need to be modeled. Accelerometers A1 and A2, which are mounted to the input 
shaft pillow block outside the gearbox (Figure 2), show a reduction of vibration amplitudes 
by at least 20 dB and do not pick up any of the input shaft modes. In addition, all four 
NASA bearing accelerometers do not pick up any vibration from this impact, not even in 
trace vibrations between 3050 H z and 4500 Hz. This suggests that the gearbox system 
under study is not coupled with the accessory drive components and permits confident 
modeling of the chosen system. 




Frequency (Hz) 


(a) 


(b) 


Figure 6: Frequency response functions for impact test to gearbox input shaft with accelerome- 
ters A8 and A9 mounted to the input shaft and accelerometers Al and A2 mounted to the input 
shaft pillow block for (a) four NASA accelerometers and (b) select OSU accelerometers. 


2.4 Impact to Rear of Gearbox Housing 

In this test an impact was made to the rear of the gearbox housing to estimate the dominant 
gearbox structural modes. Figure 3(d) shows the impact point and one of the accelerome- 
ters mounted to the gearbox (A4). The mean frequency response function from the NASA 
accelerometers is shown in Figure 7(a), and the frequency response functions from the 
relevant OSU accelerometers for this test are shown in Figure 7(b). 

Figure 7(b) shows multiple natural frequencies picked up by accelerometers mounted 
to the gearbox housing: 550 Hz, 1000 Hz, 2000 Hz, and 2800 Hz. Two of these peaks 
(2000 Hz and 2800 ITT), are evident in the NASA bearing accelerometers of Figure 7(a), 
but the two lower- frequency peaks are not apparent. This data provides an estimate of the 
primary gearbox modes and shows which of these includes significant bearing dynamics. 
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As seen before, accelerometers A1 and A2, which are mounted to the input shaft pillow 
block outside the gearbox, show a reduction of vibration amplitudes by about 20 dB and 
hardly pick up the natural frequencies at 2000 H z and 2800 Hz. This adds confidence to 
the assumption that the gearbox system under study is not significantly coupled with the 
accessory drive components, validating the modeling approach. 




Frequency (Hz) 


(a) 


(b) 


Figure 7: Frequency response functions for impact test to rear of gearbox housing with ac- 
celerometers A3 and A4 mounted to the front and rear of gearbox respectively and accelerome- 
ters Al and A2 mounted to the input shaft pillow block for (a) mean NASA accelerometers and 
(b) select OSU accelerometers. 


2.5 Impact Test Summary 

Important information was learned from these brief impulse tests. General system behav- 
ior was characterized by identifying several important natural frequencies. With thirteen 
total accelerometers, it was easy to identify the components of principal vibration at these 
natural frequencies. These tests also helped establish boundaries of the system for compu- 
tational and analytical studies. While the gears, shafts, bearings, and housing are clearly 
coupled, it is apparent that connecting input/output shafts, pillow blocks, the drive motor, 
and dynamometer structure do not contribute significantly to the dynamic behavior of the 
system. Table 1 summarizes the modes identified in this study and the most significant 
components of vibration. 


N ASA/CR— 20 1 2-2 1 7664 


16 


Table 1: Natural frequencies observed in the NASA GRC gear test rig by impact testing. 


Number 

Frequency 

Primary component 

1 

500 Hz 

Gearbox 

2 

750 Hz 

Gears/shafts 

3 

1000 Hz 

Gearbox 

4 

1800 Hz 

Shafts 

5 

2000 Hz 

Gearbox/bearings 

6 

2500 Hz 

Gears/shafts 

7 

2800 Hz 

Gearbox/bearings 

8-12 

3050 Hz -mo Hz 

Input shaft (accesso 


3 Gearbox Vibro- Acoustic Propagation Analysis 
Method Overview 

The vibro-acoustic gearbox modeling consists of four parts. They are: 1) full-fidelity finite 
element/contact mechanics model of the entire gearbox including gears, shafts, bearings, 
and the housing [18]; 2) mathematical lumped-parameter gearbox model; 3) finite element 
housing model; and 4) boundary element model of the housing [19]. 

The finite element model of the analyzed gearbox (shown in Figures 8(a) and 8(b)) in- 
cludes realistic rolling element bearings, shafts, and the housing. Other models may model 
these components analytically because of the modeling complexity. The unique contact 
algorithm seeks the gear and rolling element contact in three-dimensional space. This al- 
gorithm has been validated against experiments [25, 26, 27, 28, 29, 30, 31]. Output of this 
detailed finite element model includes dynamic bearing force, bearing stiffness, and trans- 
mission error, which are crucial information for the acoustic analysis and mathematical 
modeling. 

The analytical model of the gearbox provides fast dynamic analysis of the gearbox. 
Bearings are modeled as linear stiffnesses using 6x6 matrices computed by the finite 
element/contact analysis. Housing is included through its static condensation compliance 
matrix, computed by finite element analysis. The gear contact is modeled by a network of 
springs along lines of contact. The gear model is nonlinear and captures the non-nominal 
load distributions due to shaft tilting, misalignment, and profile/lead modification. 
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Isolated finite element housing model is established in ANSYS. The housing excitation 
source is dynamic bearing forces calculated using the finite element/contact or analytical 
model. The surface velocity of the bearing is obtained by the dynamic responses. This 
step can not be omitted because it connects the vibration analysis of the entire gearbox 
and the housing noise radiation computation in the following step. 

The full-fidelity boundary element model of the gearbox housing is established in the 
software Coustyx [32, 19]. The acoustic model employees a multipole method ([32]) to 
provide noise calculation with realistic housing dimensions. This housing surface velocity 
calculated by ANSYS is inputted into the acoustic model as its boundary condition. This 
acoustic model computes the radiated noise from the housing. 

Major steps of the gearbox vibro-acoustic analysis is depicted in Figure 9. 



(a) (b) 

Figure 8: (a) Outside and (b) inside of the gearbox at NASA Glenn Research Center. 
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Figure 9: Key steps to perform the gearbox acoustic analysis. 

4 Finite Element / Contact Mechanics Modeling and 
Analysis of the Gearbox 

4.1 Gearbox Modeling Overview 

The examined NASA GRC gearbox (shown in Figures 8(a) and 8(b)) includes a pair of 
spur gears with webbed rims, two staged shafts on which the gears are mounted, four 
rolling element bearings that support the shafts, and the housing. The gear parameters 
are listed in Table 3. Two types of rolling element bearings are used: cylindrical rolling 
bearings and deep groove ball bearings (described in Table 2). Product designations are 
SKF N205ECP and SKF 6205, respectively. The height and width of the gearbox housing 
is 279.4mm and 254.0mm with 330.2mm length. The housing wall and lid thickness is 
6.350mm. Details of the gearbox are described in [33, 24]. 

The realistic bearing model captures detailed bearing mechanics as shown in Figure 10, 
including rolling elements, races, and the cage (not shown). This detailed bearing model 
is used to determine the full 6x6 stiffness matrix between the shafts and housing. That 
matrix is used in the mathematical gearbox model to perform fast dynamic simulations, 
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Parameters (mm, degree) 


Cylindrical Roller Ball Bearing 


Number of rows 

Number of rolling elements 

Contact angle 

Pitch diameter 

Bore diameter 

Roller length 

Roller diameter 

Bearing width 

Outer diameter 

Outer diameter of inner raceway 
Inner diameter of outer raceway 
Radial clearance 
Inner race crown curvature 
Outer race crown curvature 


13 

9 

0 

0 

39.00 

38.50 

25.00 

25.00 

8.600 

7.900 

7.500 

7.900 

15.00 

15.00 

52.00 

52.00 

31.50 

34.40 

46.40 

46.30 

40.00 x 10" 3 

20.00 x 10 

o 

1 

-4 

0.520 

o 

1 

-4 

0.520 


Table 2: Single row cylindrical roller and ball bearing parameters. 


as discussed later. 

The gearbox housing is modeled by importing the full fidelity mesh established in 
commercial finite element software, PATRAN. The housing is then assembled into the 
gear/bearing/shaft system as shown in Figure 11. 

4.1.1 Contact Solver 

The contact solver of [18] seeks contact between gear teeth and bearing rolling elements 
and raceways. Mesh stiffness variation, transmission error, tooth separation, and bearing 
stiffness variation are inherently included; they are outputs rather than inputs. Transmis- 
sion error is the major excitation source in geared systems [34] so an accurate transmission 
error estimate is crucial. The finite element /contact analysis of the gearbox provides the 
required reliable transmission error estimate. 

The geometric surface descriptions of the contacting bodies must be precise to fully 
address the contact characteristics. Additionally, the contact area is narrow and travels 
over the entire body surface. Conventional finite element analysis requires a prohibitively 
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(a) (b) 


Figure 10: Cut-away finite element mesh of the radial ball bearing used in the gearbox (dimen- 
sions detailed in Table 2) and a double-row cylindrical rolling element bearing (from a helicopter 
application). 


Parameters 

Values (mm, degree) 

Number of teeth 

28 

Outer diameter 

95.25 

Root diameter 

79.73 

Facewidth 

6.350 

Module 

3.175 

Pressure angle 

20 

Center distance 

88.90 

Tooth thickness 

4.851 

Cutter edge radius 

1.270 

Linear tip relief 

0.1778 starting at 24 degrees 


Table 3: Dimensions of the spur gear pair 


refined mesh to address these problems; a complete dynamic response analysis becomes 
impossible in that case. 

The finite element/contact mechanics model used here addresses these issues by using 
a combination of the Boussinesq solution near the contacting surfaces and traditional finite 
element analysis far away from the contact zones to exploit the advantages of each. The 
details about this contact solver can be found in [35] . 

To accurately describe the contact area and pressure, the contact zone is discretized 
into many small patches (grid cells). Sufficient number of grid cells within the contact 


N ASA/CR— 20 1 2-2 1 7664 


21 



Figure 11: Assembly of the gear-bearing-shaft-housing model. 


zone is essential to obtain the correct contact pressure and load distribution. The finite 
element model of the gear pair and contact pressure on individual tooth over a mesh cycle 
are shown in Figure 12(a) and Figure 12(b). Figure 13(a) shows the finite element model 
of a double-row cylindrical rolling element bearing. Contact patches on one of the radially 
loaded cylinders are shown in Figure 13(b). 

When the gears and rolling elements rotate, the number of teeth and rolling elements 
in contact change, as do the location, size, and shape of the contact areas. These changes 
are important as they affect bearing forces, gear tooth loads, and transmission error cal- 
culations. The contact solver addresses these issues by determining and analyzing the 
instantaneous gear and bearing contact conditions at every time instant. 

This specialized finite element/contact mechanics software allows dynamic simulations 
with greater modeling fidelity than conventional finite element tools. It is validated against 
benchmark studies of complex gear dynamics problems [25, 26, 27, 28, 29, 30]. In exper- 
imental comparisons, it has proven accurate in capturing the complex tooth mesh forces 
leading to strong nonlinearity in the dynamics of single gear pairs [26], idler [29, 36], and 
planetary gears [25, 27, 28, 37]. The rolling element contact in multiple bearings has been 
validated against experiments in [31]. The shafts in the gearbox introduce system compli- 
ance, could cause misalignment, and eventually affect transmission error. The accuracy of 
the shaft models have been validated against classical beam theories as discussed later. 


N ASA/CR— 20 1 2-2 1 7664 


22 



Figure 12: (a) Mating gears in the NASA GRC gearbox; (b) contact pressure on gear teeth over 
one mesh cycle; 


4.2 Fluid Film Wave Bearing Modeling 

Bearings are critical components in geared systems. The fluid film wave bearings are 
included in the gearbox model through the stiffness and damping matrices calculated by 
the program developed by Hanford and Campbell [38, 39] . This wave bearing model uses 
a perturbation method based on the Reynold equation to calculate the stiffness, damping, 
pressure distribution, and load capacity of the fluid film. The program, however, is limited 
to be two-dimensional by excluding the bearing tilting motion. 


4.3 Rolling Element Bearing Modeling and Analysis 

Theoretical bearing models [1, 40, 41, 42] make different assumptions to formulate the load- 
deflection relation. These assumptions include different race elasticity, different property of 
the rolling element contact, and ignoring microgeometry dimensions of bearing raceways. 
Figure 14 shows the nonlinear stiffness-load relations of the cylindrical and radial ball 
bearings calculated using the Harris [1], Gargiulo [2], and While [3] models. Significant 
discrepancy is present among them. 

The method developed to determine bearing stiffness does not make any assumptions 
about the load-deflection relation. Instead, it calculates the 6x6 bearing stiffness matrix 
by partial derivatives of applied forces and moments related to six degrees of freedom. This 
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Figure 13: (a) Finite element model of a double row cylindrical bearing (outer race is removed); 
(b) contact pattern on one of the loaded cylinders. 


stiffness determination approach is detailed in [31]. 

Radial and axial stiffness of radial ball bearing in [4] calculated by the proposed method 
under axial preloads is compared against the experiment conducted by Kraus et. al. in [4] 
as shown in Figure 15. Because of the lack of information on internal radial clearance, two 
clearance values of 0 and 0.01 mm are used to compute stiffness. Radial and axial stiffness of 
self-aligning ball bearing in [5] calculated by the proposed method is compared against the 
Royston and Basdogan [5] experiment under radial and axial preloads, respectively (Figure 
16). Good agreement is evident between the proposed method and the measurements for 
both bearings. 

The number of rolling elements in contact oscillates when they rotate, which indicates 
bearing stiffness is periodic about the ball pass frequency. The ball pass period T = 
~ jj t is the amount of time between one rolling element leaving a reference point and 
the next rolling element arriving. Here, z is the number of rolling elements, u denotes the 
shaft speed, and quantities D e and D p denote the diameter of rolling elements and pitch 
diameter of the bearing, respectively. 

Figures 17(a) and 18(a) show radial stiffness of the cylindrical and radial ball bearings 
in a ball pass period with 1000 N applied load. Radial stiffness of the cylindrical bearing 
fluctuates with 16% deviation from the mean stiffness, while the number of rolling elements 
in contact alternates between four and five over the ball pass period. Figures 17(b) and 
18(b) show tilting stiffness of the cylindrical and ball bearings when 1 Nm moment is 
applied. These stiffness fluctuations can excite gearbox vibration. 


N ASA/CR— 20 1 2-2 1 7664 


24 



600 



Figure 14: Radial stiffness of examined cylindrical bearing and radial ball bearing vs. applied 

radial loads calculated by the Harris (— ) [1], Gargiulo (• — •— ) [2], and While ( ) [3] models. 

The While [3] model is modified to use ^ to calculate the stiffness instead of 


4.3.1 Cross-Coupling Bearing Stiffnesses 

Traditionally, diagonal stiffness matrices are used to represent rolling element bearings. 
These stiffness matrices, however, include both diagonal and off-diagonal terms. Equation 
1 demonstrates the stiffness matrix structure. The quantities k xx , k yy denote radial stiff- 
ness. k zz denotes axial stiffness. The quantities kg x g x , kg y g y denote tilting stiffness, which 
prevents the tilting motion of the shafts. The off-diagonal stiffness falls into four cate- 
gories: the coupling between radial and rotational displacements (k x g x , k x g y , k y e x , k y g y ), 
the coupling between radial and axial displacements (k xz , k yz ) , the coupling between axial 
and rotational displacements ( k z g x , k z g y ), and other coupling terms ( k xy , kg x g y ). Rolling 
elements are free to rotate in 8 Z direction so that stiffness in 9 Z direction and other related 
matrix components are zeros. The stiffness matrix is symmetric because rolling element 
bearings are conservative systems. 

k'f 


K = 


k X y 

k X z 

k;r() :l: 

k X 6y 

0 

kyy 

h 

l^yz 

k y0 x 

k y6y 

0 


kzz 

k z e x 

k Z 9y 

0 



k e x e x 

k e x o y 

0 

Symmetric 



ke y e v 

0 


( 1 ) 


Cross-coupling terms in the stiffness matrix indicate interactions between radial, axial, 
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Figure 15: Comparison between the proposed method with zero (•••□•••) and 0.01 mm (• — 
□ • — ) radial clearances and Kraus et al.’s [4] experiment (— o — ) for radial and axial stiffness 
of the ball bearing in [4] under axial preloads. 

and tilting motions of rolling element bearings. They demonstrate the coupling between 
the shaft tilting motion, the flexural motion of the structure connected to the outer race, 
and the shaft radial and axial motions. The effects of cross-coupling terms on the gearbox 
vibration transmissibility through rolling element bearings are investigated. As the primary 
excitation source, gear transmission error is an important measure of gearbox vibration. 
Figure 19 shows the spectra of dynamic transmission error in the frequency range from 1500 
to 4000 Hz from numerical torque impulse cases. Bearing models with fully-populated and 
diagonal stiffness matrices are compared as shown in Figure 19. Differences in resonant 
frequencies and amplitudes are evident between these two bearing models. This stresses 
the significance of the cross-coupling stiffnesses. 

4.4 Gear Transmission Error 

Transmission error is the major excitation source in geared systems. Accurate transmis- 
sion error estimate is crucial. With the precise contact solver, the finite element/contact 
analysis of the gearbox provides reliable transmission error estimate. 

Transmission error is computed according to the tooth mesh deflection TE = {x\ — 

X 2 )sin(a) + (y\ — y 2 )cos(a ) + r\6\ + r 2 & 2 , where Xj, jy, 6i, i = 1, 2 are the coordinates of 
the mating gears as shown in Figure 20. a is the pressure angle and r % , * = 1,2 are the 
radii of the meshing gears. This formulation includes the shaft and bearing compliance. 
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Figure 16: Comparison between the proposed method (•••□•••) and Royston and Basdogan [5] 
experiment (— o — ) for radial and axial stiffness of self-aligning ball bearing in [5] under radial 
and axial preloads, respectively. 


Transmission error of this gear pair without shaft or bearing compliance has been com- 
pared among Program X, Load Distribution Program, NASA DANST [43], and the current 
approach in Figure 21. Program X is multi-body dynamics software that is used by in- 
dustries worldwide. We are not free to state its name because of license restrictions for 
academic use. The agreement on the peak to peak amplitude is reasonable. Differences 
among the mean amplitudes are present. These differences are mainly caused by different 
rim models these programs have used. Gear rims introduce compliance into the system, 
leading to high amplitude of the mean transmission error. The current model includes the 
realistic rim (as shown in Figure 12(a)). Others model the rims differently by excluding 
the rim shoulders. Thus, their estimates of transmission error are lower. The Harris map 
of transmission error at various torques is shown in Figure 22. The minimum transmis- 
sion error without including the shaft compliance is at 67.79 Nm torque, which matches 
the torque the gear teeth are modified at. This further validates the transmission error 
estimate. 

In addition, the minimum of the peak to peak value of transmission error is at lower 
torque (57.62 Nm) when flexible shafts are included as shown in Figure 22. This suggests 
shaft compliance needs to be considered to estimate transmission error and modify gear 
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(a) (b) 

Figure 17: (a) Radial and (b) tilting stiffness of the cylindrical bearing over a ball pass period. 
The bottom figures show the number of rolling elements in contact over a ball pass period. The 
applied load and moment are 1000 IV and INm, respectively. 

teeth. 

4.5 Shaft Modeling and Validation 

The long shafts in the gearbox introduce system compliance, could cause misalignment, 
and eventually affect transmission error. The accuracy of the shaft models are important. 

The shaft bending and torsional deformations are compared against classical beam 
theories. The shaft bending model is considered as the elastic beam with a concentrated 
load applied at the gear location. The beam bending boundary conditions are chosen as 
simply-supported at each end of the shaft. The transverse deflection y along the shaft is 
calculated as 

Pb[x 3 — (L 2 — b 2 )x ] 3 / 2 

V — qeTl ’ x < a 

(2) 

r _ < R iut~RL) 

1 ~ i 

where the parameter P is the applied concentrated force. The parameters a, b are the 
distances between one shaft end and the location where P is applied. The quantities 
E, /, L denote the Young’s modulus, moment of inertia, and shaft length. The quantities 
Rout , Rin denote the shaft outer and inner radii. 

The analytical model to calculate the shaft torsional deflection has the clamped-free 
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(a) (b) 


Figure 18: (a) Radial and (b) tilting stiffness of the ball bearing over a ball pass period. The 
bottom figures show the number of rolling elements in contact over a ball pass period. The 
applied load and moment are 1000 IV and INm, respectively. 

boundary condition. A torque T is applied at the free end of the shaft. The torsional 
deflection 0 is calculated as 

9 = Tx/GJ 

o = P) 

J = § - K.) 

where the quantities v, J denote the Possion’s ratio and the second moment of inertia. 

Figures 23(a) and 23(b) show the shaft bending and torsional deformations calculated 
by analytical solutions (solid line) and finite element results (square marker) at various 
torques/forces. The finite element results of shaft bending and torsional deformations 
agree with the analytical predictions. 

In addition, as shown in Figure 4.5, the bending deformation is the same order of 
magnitude of transmission error near the operating torque. Shaft deformation is significant. 

Thus, including shafts would increase the overall accuracy of gearbox modeling. 
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Figure 19: Numerical torque impulse response of gear dynamic transmission error with fully- 

populated (— ) and diagonal ( ) stiffness matrices of the rolling element bearings mounted in 

the examined gearbox based on [6]. The input torque equals 84.74 ATn. 


y, 


— v, 

e 2 

Figure 20: Coordinates of the examined gear pair. 

4.6 Dynamic Analysis and Correlation with Experiments 

Modal analysis is performed numerically on the gearbox model by applying a torque im- 
pulse at the input shaft and measuring the dynamic response. Numerical impulse test 
provides natural frequencies and mode shapes of the gearbox. Results of the impulse test 
are correlated with experimental measurements conducted by NASA personnel [44] and 
the Ohio State research group. 

The frequency spectrum of dynamic transmission error and shaft displacement for the 
speed range from 0 to 8000 Hz are shown in Figures 25(a) and 25(b). The computed 
natural frequencies agree with experiments as compared in Table 4. NASA accelerometers 
identify natural frequencies near 3000 to 4500 Hz, but the coherence is relatively low at 
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Figure 21: Static transmission error of the gear pair without the shaft and bearing compliance. 
The results are calculated by finite element (— ), Program X (— • — ), NASA DANST (•••), and 
Load Distribution Program ( ). The torque equals 79.09 Nm. 

these higher frequencies. Measurements show several modes near 6500 to 7500 Hz that 
are difficult to resolve, while simulations predict only one natural frequency at 6856 Hz. 

The mode shapes of these natural frequencies include mesh deflection modes, shaft modes, 
housing modes, and coupled modes. 


Finite Element 

Experiment 

759 

750 

1942 

2000 

2195, 2283 

2500 

3803, 4449, 4654 

3000-4500 

6856 

6500-7500 


Table 4: Natural frequencies predicted by numerical impulse tests and measurements (Hz) 
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Figure 22: Peak to peak amplitude of static transmission error of the gear pair with (— x 
and without (— o — ) shafts at different torques. 

5 Analytical Gearbox Dynamic Modeling and Anal- 
ysis 

Mathematical modeling of the gearbox provides fast dynamic analysis of the gearbox. 
Dynamic analysis, especially at multiple operating speeds, requires much computational 
effort using the finite element model. Analytical model gives the natural frequencies and 
mode shapes, and dynamic response at various operating speeds. 

This analytical model employees crucial information on transmission error, mesh stiff- 
ness, bearing stiffness, and housing compliance from the finite element analysis. 

5.1 The Gear Pair 

The gear pair model is constructed following the literature [45, 46, 47, 48, 49, 50, 51], in 
which detail discussions into the validity of the modeling can be found. The model consists 
of two gears mounted on shafts. Each gear body is combined with its supporting shaft 
into a single rigid body. These gear-shaft bodies are each mounted on up two bearings 
placed at arbitrary axial locations. Figure 26 shows the gear model and the bases. A fixed, 
right-handed, orthonormal basis {E} = {Ei, E2, E3} is oriented such that Ei is parallel to 
the line of action of the gear mesh. The origin is on the rotation axis of the pinion body, 
midway in the active facewidth. This dimensioning allows arbitrary axial positioning of 
meshing gears with different facewidth because any inactive facewidth is treated as part of 
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(a) (b) 


Figure 23: (a) Input shaft bending deformation calculated by analytical beam theory (solid line) 
and finite element method (square marker) with simply-supported boundary conditions under 
various input torques; (b) Input shaft torsional deformation calculated by analytical beam theory 
(solid line) and finite element method (square marker) with clamped-free boundary conditions 
under various input torques. The shaft has uniform outer diameter (30.23 mm). 


the shaft. The translational ( x p ,y p ,z p ) and angular (</> p , 9 P , (3 p ) coordinates of the pinion 
body are assigned to translations along and rotations about Ei, E 2 , and E 3 , respectively. 
The translational and angular coordinates of the gear body follow similarly with subscript 
g. Body-fixed bases {e p } = {e^, e^, e|} and {e 9 } = {e^ 1 , e®, e^} for the pinion and gear are 
adopted. Positive axial quantities are measured along E 3 from the dashed line in Figure 
26. 

The pinion translational and angular velocity vectors are 


Fp — XpEi -\- ijpEi2 T Z P E 3 , 


LJp — 





/3p + 



(4) 


where Q p is the constant angular rotational speed of the pinion. The velocity vectors for 
the gear are identical except with components for the gear. 

The pinion body is supported by two bearings at points A p and B p . The axial positions 
of these bearings measured along E 3 are L p and L ^ . The pinion bearing deflection vectors 
at point A p and B p are the relative deflections of points A p and B p with respect to ground, 
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Figure 24: Shaft bending deformation under various input torques. 


giving 

dp = | dp [Lp epj + Xpj Ei + | (f>p (e.p — L p ) + y p j E2 + z p E 3 , 

(5) 

d p = \0 p ( Lp — e p ) + x p ] Ei + [</> p (e p — + y p ] E2 + £ p E 3 . 

The bearing deflections for the gear follow similarly. The bearings resist tilting as well. 
The angular deflection of the pinion body bearing at A p is 

7 p = 1 + 0 p E 2 + ^ p E 3 . (6) 


The angular bearing deflection at point B p is identical to Eq. (6) for rigid shafts. The 
bearings are isotropic in the Ei — E 2 plane. At point A p , the bearing stiffness matrix K p 
as given in Eq. (1) is fully-populated, where the equality of stiffness in the two translation 
directions is evident. The bearing translational and angular displacements combined are 


r = 



7 


Similar definitions follow for point B p and for the gear body. 


The gear mesh interface is modeled by a series of springs along the nominal lines of 
contact for no mesh deflection. These lines change as the gears rotate. Each spring acts 
at a point denoted by Ci. When the gear bodies deflect, the contact points on the pinion 
separate or compress against the contact points on the gear. The difference between 
the position vectors of the contact points on the pinion and gear gives the relative mesh 
deflection vector at Ci. The projection of the relative mesh deflection vector on the tooth 
surface normal gives the relative compressive deflection at the zth contact point. The 
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Figure 25: Numerical impulse test results of (a) dynamic transmission error and (b) the input 
shaft horizontal displacement of the gear-bearing-housing system within speed range from 0 Hz 
to 7000 Hz. The applied torque is 79.09 Nm. 


relative compressive deflection is 

5i(q, t) = \ [e p - Ci{t )] 9 P + [ a(t ) - e g \ 6 g - x p + x g + + /3 p r p + /3 g r g l cos ip 

r ' ' ' “ i ( 7 ) 

- | [k{t) + hi] 0 P + [(r p + r g ) 0 g tan<f> - bi(t)] + z p - z g + cp p r p + (p g r g j sin^, 

where r p and r g are the base radii, 4* is the transverse operating pressure angle, and ip is 
the base helix angle. The vector q comprises generalized coordinates 

q = ( (j>pi fipi x pi Vp-i z pi 0 91 figi x gi Vgi z g ) (8) 

pinion gear 

The axial position of a contact point is Cj(i) measured from the origin along E 3 , and the 
transverse position of a contact point is bi(t) measured from the origin along — Ep. They 
are known functions of time determined by the contact line progressions as the gears rotate. 
Micron- level deviations of the tooth surface from an involute at any contact point i, such 
as from gear tooth surface modifications and manufacturing errors, are denoted by hi. 
Figures 26 and 27(b) depict these quantities. 
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Figure 26: Analytical model of the gear pair. The parameters are defined in [7]. The dashed 
line is at the center of the active facewidth. 


The kinetic and potential energies are 


T = ^ (cjp J p u p + mjj g Ug + TpUipip + Y T g m g r g ) , 

1 r iV T K®iy + r® T K,®rJ + + rJ T *'®r B 


v= 2l 


■ G 1 - p 


-G*-p 


(9) 


+ 


n(t) 

y 

2 ^ 




i — 1 

where ki(q,t) is the *th contact stiffness, and n{t) is the number of contact segments at 
an instant t. These quantities change as the gears rotate, hence the time dependence. 
The inertia tensor of the axisymmetric pinion body is J p = diag [Jp, Jp, Jp\ with similar 
definition for the axisymmetric gear body. 

Lagrange’s equations of motion for unconstrained generalized coordinates follow after 
substitution of equations Eqs. (4) through (7) into the energy expressions Eq. (9). In 
matrix form they are 


M s q s + Dq s + fl p G s q s + [K s (q, t) - OpC s ] q = F(q, t) s , (10) 

The vector F includes external loading; the driving and absorbing torques and tooth surface 
modifications hi appear here. The matrix K represents the system elasticity with losses 
contained in the modal damping matrix D. Tooth surface modifications hi are neglected 
in K because hi <C bi(t), and the hi appear as additions to bi(t). The terms that arise from 


N ASA/CR— 20 1 2-2 1 7664 


36 



Figure 27: (a) Distributed spring network over a contact line with the local and bulk stiffnesses, 
k c l(t ) and kb. (b) Local k c l(t) and bulk kb stiffnesses are combined into contact stiffness fcj(q, t) 
by Eq. (14). 

the constant rotation speed are contained in the gyroscopic matrix G and the centripetal 
acceleration matrix C. Individual elements of M, K, G, C, and f are given in [7]. 

Following [52, 50] the nominal contact lines are discretized into n(t) segments of equal 
length l(t), as shown in Figure 27(a). Each contact point Ci is positioned at the center of 
its segment. As the contact lines progress with gear rotation, the total number of segments 
n(t) and the length of a segment l(t) change. Each contact line has a specified number 
of segments. This discretization is based on the nominal lines of contact with no gear 
deflections. 

Each contact spring is attached to its contact point Ci. The stiffness fcj(q, t) of contact 
springs are obtained by considering two separate categories of tooth deflection: local (ej) 
and bulk (6b). Discussion of this categorization can be found in [52, 49, 53]. The local 
deflection represents the Hertz contact deflections. The associated local stiffness is k c l(t), 
where the constant k c is the local stiffness per unit contact length. The bulk deflection 
represents all deflections except local deflection, and those include gear blank deflection, 
tooth bending, shear, etc. Because the Hertz contact deflections are localized and far 
enough from the bulk deflections, the bulk deflection is assumed to be the same for all 
contact segments. The bulk stiffness kb is assumed constant. The bulk spring is in series 
with the local springs, so the total deflection at the zth contact point Ci is 

$i = e i + $b (11) 

The mesh force F equals the sum of all forces carried by the local springs and also the 
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force carried by the bulk spring due to the series connection. The mesh force is 


n n 

F = Y J F l = k c l(t ) ^ e t H(ei) = k b 5 b , (12) 

i=l i = 1 


H(ti) 


l; ei>0 
0 ; e» < 0 


(13) 


is the Heaviside function that represents the contact condition at each contact spring. Use 
of Eqs. (11) and (12) reduce the network of local and bulk springs into n(t) contact springs 
(hi, i = 1,2, ... ,n) in parallel, as shown in Figure 27(b). The zth contact stiffness is given 
by 


fci(q,f) 


k b k c l(t)H(ei) 

n(t) 


(14) 


k b + k c l{t) y^ H(ej) 


i = 1 

One can then calculate the zth potential energy ki5f / 2 stored in the gear mesh in Eq. (9) 
from the contact force F t = kjS t . 

The local stiffness per unit length k c and the bulk stiffness k b are parameters of the 
gear pair determined by the contact mechanics and elasticity of the gears. These constants 
can be approximated analytically [52] or semi-analytically [50] by assigning certain types 
of stiffnesses such as Hertz contact, tooth bending, shear, etc. to k c and k b . A different 
approach is used in this work, where k c and k b are solved for from the deflections obtained 
from an external analysis tool for computational static analyses. In this case, k c and k b are 
numerical values that best fit the deflection obtained from finite element contact analysis 
of gears. 

The following stipulations simplify the algebra to find k c and k b . 1) The tooth surface 
is perfectly involute, that is, hi = 0 for all i; 2) All degrees of freedom are constrained 
to be zero except the pinion rotation /3 P ; and 3) A specified moment about E 3 is applied 
to the pinion. With these stipulations, the deflections at all contact points are identical, 
that is, 8i = 5 for all i in Eq. (7). Consequently, all points are in contact; Hi = 1 for 
all i. The subscript i of ki is unnecessary because when all segments are in contact, 
h = k 2 = ■ ■ ■ = k n . Use of static equilibrium, Eq. (11), and Eq. (14) gives 


k b k c L(t)\ ’ 

where S(t) is the static transmission error, L{t) = n(t)l{t ) is the total contact line length 
at an instant t, and F is the constant mesh force obtained from the known applied torque. 
The two unknowns {k c , k b ) are solved using the data from finite element analysis results 
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at two instances {(f(fi), L(t\)} and {Sfo), Lfo)} within a mesh cycle. To increase accu- 
racy, values at these two instances are calculated from averages of the four points where 
transmission error is highest (giving the values for the first instance) and the four points 
where transmission error is lowest (giving the values for the second instance). 

The gear mesh model is validated by finite element analysis. Figure 28 shows trans- 
mission error comparison between the finite element and the analytical model for a helical 
and a spur gear pair. 



(a) (b) 

Figure 28: Static transmission error from the analytical (solid line) and finite element (circles) 
model, (a) A helical gear pair. Quadratic tip relief starting at a = 28 deg and root relief at 
a = 27 deg. Tip relief, root relief, and circular lead crown are 10 /zm. The applied torque is 
200 N-m. (b) Spur gear pair in [8]. Linear tip relief starting at a = 23.6 deg with amplitude 10 
jim. Circular lead crown is 5 /zm. The applied torque is 340 N-m. 

5.2 Incorporating Housing Compliance with the Analytical 
Gear/Shaft Model 

The equations of motion for the housing is given by 

Mi/q# + K/f-q# = F# (16) 

where M# and K h are the mass and the stiffness matrices with degrees of freedom q u 
and the force vector f#. Figure 29 helps explain the approach. The gear/shaft equations 
of motion are expanded by the addition of bearing degrees of freedom q;, and the housing 
equation of motion by Partitioning of (16) to separate the degrees of freedom where the 
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bearings are attached from the rest of the housing giving 


" M c 

0 

( 4 ) + " 

0 

M h 

V <lh ) [ 


K r 


Kh 


K ch 

K h 


The new equations of motion are given by 


" M, 

0 

42 : 

Co 

+ 

1 

0 

M b 

\ / 


K, 


*1 


K sb 

K, 


q c 
q h 


q s 

q& 



(17) 


(18) 


where, the matrix K s ;, includes the coupling between the bearing degrees of freedom and 
the system. M c and K c are the mass and stiffness matrices, q c is the degrees of freedom, 
and f c is the force at the connection points. The bearing ends are massless, M& = 0. 

The rest of the housing mass, stiffness, degrees of freedom and forcing are contained in 
Mfr, K/i, q^, and f^. The connection between the bearings and the housing requires that 
at the interface displacements are equal, q& = q c , and the forces are transmitted f b = — f c . 
Eliminating q c using these connection requirements gives 


" M s 

0 

0 


^ q s ^ 


" K s 

K. s6 

0 


^ q s ^ 


/ 

is 

\ 

0 

m 6 

0 


q b 

+ 

K b 

k c + k 6 

K ch 


q b 

= 


0 


0 

0 

M h _ 


V ^ J 


0 


K h _ 


V / 


V 

f h 

/ 


This is the general form that connects the gears, shafts, and bearings to the full housing 
model which includes the flexibility as well as the inertia. If the housing inertia is neglected, 
i.e. static condensation, M/, = M 5 = 0 and there are no external forces on the housing 
f h = 0, the last row of (19) becomes q/ t = — K^K^qf,. Substitution of q^ into the second 
row of (19) gives q& = — [K& + K C -K*K,X] K^q,. Substitution of q& into the 
first row of (19) gives 


-K ch K^K^ h } 

s 

V 

K/ 

This equation incorporates the housing flexibility, without added degrees of freedom and 
the housing inertia, into the equations of motion of the gears, shafts, and bearings. To this 
end, the calculation of K c , K c ^, and K bl separately is unnecessary. A shorter path can be 
taken by evaluating the housing stiffness Kf at the connection interface using influence 
coefficients of the housing (16). This is discussed in the following subsection. 


}q s = 0 (20) 


Msqs + — [K b + K c 
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Figure 29: Description of the connection between the gears, shafts and bearings to the housing 
model. 

5.2.1 Including the Housing by Using Influence Coefficients 

The influence coefficient matrix C f can be solved from (16) Cjf c = q c Solving for 
from (17) gives 

Cj 1 =K f = K c ~K ch K- 1 K^ h (21) 

which allows use of the influence coefficient matrix in (20). 

The procedure to obtain elements of CJ xn consists of applying a unit force to each 
degree of freedom at the connection points f c one at a time, and measuring the deflection 
of all the degrees of freedom at the connection points q c . Explicitly, elements in the i th 
column, Cj i><J> , are given by the displacements q^ J> , calculated by (16) when there is 
a unit force at = 1, where j = 1 , ,n. 

5.3 Dynamic Analysis Using the Analytical Model 

Table 5 shows the natural frequencies of four systems: 1) cylindrical/bail bearings with 
housing flexibility, 2) wave bearings with housing flexibility, 3) cylindrical/ball bearings 
without housing flexibility, and 4) wave bearings without housing flexibility. Figure 30 
shows the mesh deflection mode from the analysis with cylindrical/ball bearings with 
housing flexibility. The mesh deflection mode is an important mode because the tooth 
pass excitation excites primarily the mesh deflection mode. This mesh deflection mode is 
demonstrated in experiments, which is strongly excited by the 2 nd harmonic of tooth pass 
frequency at 4200rpm. This mode shows up strongly in our vibration speed maps. It is 
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highlighted in Table 5. The mesh deflection mode persists when bearing type is changed, 
and when the housing flexibility is included. 

Linear dynamic analysis is performed using the mean system stiffness matrix obtained 
by averaging the quasi-static time-dependent system stiffness over a mesh period. The 
harmonic excitation is approximated from the first ten harmonics of the quasi-static dis- 
placement vector. Figure 31 shows the dynamic peak-to-peak transmission error from 
the four systems mentioned above. All systems have similar dynamic transmission error. 
Figure 32 shows the dynamic bearing reactions of the four systems. Wave bearings seem 
to produce higher dynamic bearing loads. Use of the housing flexibility moves the natu- 
ral frequencies and thus the resonant peaks, but does not alter the peak amplitudes in a 
significant manner. 


Cylindrical Roller Ball Bearing 


With Housing 

Without Housing 

With Housing 

Without Housing 

576.6 

975.5 

875.2 

2279 

954.6 

1000 

1988 

2535 

987.7 

1597 

2246 

3173 

1515 

1711 

2387 

3209 

1599 

1868 

2749 

3480 

1793 

1977 

3003 

3549 

1915 

2736 

3150 

3982 

2075 

2848 

3373 

4383 

2356 

3910 

3514 

4473 

2722 

- 

3957 

- 

3893 

- 

4392 

- 


Table 5: Natural frequencies of the NASA gearbox with cylindrical roller/ball bearings, wave 
bearings, with housing flexibility and without housing flexibility. The mesh deflection mode 
natural frequencies are highlighted. 


N ASA/CR— 20 1 2-2 1 7664 


42 



Figure 30: Graphical representation of the 12th mode at 3893 Hz (mesh deflection mode) from 
the system with ball/cylindrical bearings with housing. 

6 Acoustic Radiation Gearbox Modeling and Anal- 
ysis 

The boundary element model of the gearbox housing is shown in Figure 33. This model 
includes 50,401 tetrahedron elements. Each element has nearly uniform length in order 
to have good aspect ratio. This model uses the fast multipole method (detailed in [32]) 
to calculate the gearbox sound radiation. This method provides fast acoustic calculation 
with a realistic system modeling. Given the complexity of housings, the great increase in 
the number of degrees of freedom admits solutions for realistic models not possible with 
conventional boundary element software. This method is implemented in Coustyx [19]. 

6.1 Model Validation and Mesh Convergence Study 

This boundary element model of the gearbox has been validated against theoretical solu- 
tions of the sound field with a monopole. A monopole is a source which radiates sound 
equally well in all directions. The simplest example of a monopole source would be a 
sphere whose radius alternately expands and contracts sinusoidally. The far-field pressure 
radiated by a monopole may be written as 

P(r,<?,t) = ! r^e i ( wt - fcr ) (22) 

47rr 

where p is the fluid density, c is the speed of sound, k is the wave number, oj is the 
frequency, and r is the local distance from source. Q is the complex source strength, which 
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Figure 31: Peak-to-peak dynamic transmission error of four systems. Analysis with roller ele- 
ment bearings are marked by B, analysis with wave bearings are marked by W, analysis including 
the housing flexibility is marked by H, analysis without the housing is marked by nH. 

is a constant number defined in [54], 

The sound field excited by a monopole is not affected by the gearbox geometry, which 
allows simulation results to be compared against the theoretical solution of a monopole 
for an arbitrary geometry. In other words, the noise radiation from the gearbox surface 
equals to the theoretical solution of a monopole in the free space (shown in Figure 34). 

The velocity field caused by the point source (monopole) in free space at gearbox surface is 
calculated theoretically and used as the boundary condition for the actual housing model. 

Optimal element length has been carefully chosen through mesh convergence study. 

The relative error of the computed pressure compared to the theoretical value in (22) 
with various element length is shown in Figure 35. Result accuracy is improved when 
element length is small. However, when element length is smaller than 8.636 mm, the 
result accuracy is disrupted as round off error becomes significant. The optimal element 
length is between 8.890 mm and 10.16 mm with less than 2% relative error. The mesh 
convergence study is important as the element length affects the solution accuracy. 

Sound pressure of the radiated noise calculated by Coustyx and the theoretical solution 
is compared in Figure 36 for the speed range from 500 to 3500 H z using 9.000 mm element 
length. Excellent agreement is present. 
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Gear bearing A 



Gear bearing B 



Pinion bearing A 



Frequency [Hz] 


Pinion bearing B 



Frequency [Hz] 


Figure 32: Dynamic bearing forces. Analysis with roller element bearings are marked by B, 
analysis with wave bearings are marked by W, analysis including the housing flexibility is marked 
by H, analysis without the housing is marked by nH. 


6.2 Sound Pressure and Power Computation Using Transfer 
Functions 

The computation of sound pressure and power of the radiated gearbox noise uses a network 
of transfer functions from unit bearing loads to sound pressure. Transfer functions describe 
the noise radiation property of the gearbox housing. The computation of these transfer 
functions is independent of the vibration analysis of the bearing reaction forces. Thus, 
these transfer functions are applicable for any type of bearings or gears. 

The procedure to calculate transfer functions consists of applying a unit dynamic force 
to each degree of freedom ( x,y,z,9 x ,6 y ,0 z ) at every bearing one at a time in the finite 
element model, measuring the velocity at the gearbox surface for various speeds in the 
finite element housing model, and computing the sound pressure of the radiated noise at 
specified locations in the gearbox environment using the computed velocity as the boundary 
condition in the boundary element model. Figure 37 shows a few calculated transfer 
functions at four different bearings. Sound pressure at arbitrary locations outside the 
gearbox is determined by the superposition of these transfer functions and multiplying the 
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Figure 33: Boundary element model of the gearbox established in Coustyx. 

fluctuating amplitudes of dynamic bearing forces at each degree of freedom of individual 
bearings. 

Sound power is a better measure of the gearbox noise than sound pressure as it is inde- 
pendent of the geometric locations where noise is measured. Sound power measurements 
in experiments, however, require much effort. This study calculates the sound power using 
the aforementioned transfer functions. 

To obtain accurate estimate of sound power, many numerical microphones are outside 
the gearbox. Those microphones are mounted on a measurement surface enveloping the 
gearbox, in this study, a sphere. According to the international standard ISO 3745, 20 
microphones are used. The sound power is computed by integrating the sound pressure at 
these “microphones” in the sphere. 

The sound power level, L w is computed as: 

§ 

L w = L p + 10log w (—)(m dB) + c\ + c 2 
so 

« = dB) (23) 

c 2 = dB) 

where s is the area of the measurement surface and so = 1 m 2 . B is the barometric air 
pressure during measurements, in Pascals; Bq is the reference barometric pressure. 9 is 
the air temperature during measurement. L p is the weighed surface sound pressure level 
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Theoretical Boundary Condition Excited by Monopole 


Figure 34: Theoretical solutions used to validate the boundary element housing model. Noise 
radiated from the housing with monopole velocity field at the gearbox surface as the boundary 
condition equals to that with monopole in the free space. 


over the measurement surface in decibels. 


L p = 10 log w {(±[J2 10 ai(i - +W) ])}in dB 


L pi = 20log w ( — ) 
Po 


(24) 


where N is the number of microphone positions. W is the weighting function applied 
by the filter at the frequency of analysis, pi is the root mean square pressure in Pa and 
Po = 2 x 10“ 5 Pa. 


6.3 Radiated Noise Correlation with Measurements 

The calculated noise is compared against the measurements taken at NASA GRC gearbox. 
Two microphones are mounted to measure sound pressure in experiments at 790 nun 
directly above the gearbox, and are separated horizontally by 180 mm as shown in Figure 
38. 

Here, only the noise transmitted through rolling element bearings are compared against 
measurements because experimental data on sound pressure with wave bearings is not 
available. Figure 39 shows the time averages of the mean squared sound pressure ( p)\ v = 
Sb ( Pb)av a t microphone 1 estimated by measurements and computed, where ( Pb)av 2 is 
the addictive measure of sound pressure associated with the frequency band b from 0 to 
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Boundary Element Length (mm) 


Figure 35: Effects of boundary element length on the relative error of calculated sound pressure 
compared to the theoretical solution. 


8000 Hz. For each gear mesh frequency, two different experimental results are shown: 
the minimum sound pressure including only mesh frequency harmonics; and the maxi- 
mum sound pressure including mesh frequency harmonics and other frequency components. 
These frequencies are the side bands caused by the shaft rotation frequency, belt splice 
pass frequency, etc. The simulated sound pressure is within the experimental data range. 

Figure 40 shows the frequency spectrum of the measured and simulated sound pressure 
when mesh frequency equals 2000 Hz. Good agreement is evident between the measure- 
ments and simulations at mesh frequency harmonics. Sidebands around the mesh frequency 
harmonics are present in the measured data. These sidebands are not included in the com- 
putational model because the causes are not known. Including these sidebands in the 
vibration/acoustics model would provide better agreement with the measurements. 
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Frequency (kHz) 


Figure 36: Sound pressure at the NASA microphone 1 location calculated by Coustyx (□) and 
theoretical models (— ). 



Figure 37: Sound pressure transfer functions when unit dynamic loads are applied at bearings. 
Six transfer functions are generated per each bearing along x, y, z, 9 X , 9 y , 9 Z directions. The input 
torque is 79.09 Nm. 
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Figure 38: Microphones (1 and 2) mounted above the NASA GRC gearbox. 



Figure 39: The time average of the experimental ( ) and calculated (— ) mean squared sound 

pressure at the microphone 1 location within mesh frequency range from 500 to 3000 Hz. The 
applied torque is 79.09 Nm. 
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Figure 40: Frequency spectrum of the measured (left) and simulated (right) sound pressure at 
the microphone 1 location when mesh frequency is 2000 Hz. The applied torque is 79.09 Nm. 


6.4 Noise Radiation Properties with Different Bearings 

This section investigates the effectiveness of fluid film wave bearings in reducing noise. To 
that end, sound pressure and power of the radiated noise are compared between fluid film 
wave and rolling element bearings. 

Sound pressure at the microphone 1 location excited by the first (left) and second 
(right) mesh frequency harmonics is compared between rolling element and wave bearings 
in Figure 41. The sound pressure of wave bearings is generally higher than with rolling 
element bearings. The difference is frequency dependent. When the second mesh frequency 
harmonic is near 3803 Hz (dominant mesh deflection mode), sound pressure with rolling 
element bearings is similar in shape and amplitude to that with wave bearings. The noise 
radiation is independent of bearing types at this frequency. 

The sound power radiated from the gearbox is also compared between wave and rolling 
element bearings in Figure 42. The sound power level of the gearbox with rolling element 
and wave bearings is also speed-dependent. The radiated gearbox sound power with rolling 
element bearings is lower than that with fluid film wave bearings in general within the speed 
range of interest. There is no obvious benefit to replace rolling element bearings with wave 
bearings to reduce the gearbox noise within the speed range of interest, given the two 
dimensional wave bearing model used [38, 39]. 
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Mesh Freq (Hz) Mesh Freq (Hz) 

Figure 41: The frequency spectrum of the measured (left) and computed (right) sound pressure 
at the microphone 1 location when mesh frequency is between 1000 Hz and 2500 Hz. The 
applied torque is 79.09 Nm. 

7 Summary and Conclusions 

This work investigates noise radiation from a gearbox with different types of bearings. A 
multi-dynamics model of the examined gearbox is established to address the vibro-acoustic 
propagation of the gear dynamics. It includes: 1). full fidelity finite element gearbox 
model; 2). lumped-parameter gearbox model; 3) finite element housing model; and 4). 
boundary element housing model. The procedure of the noise radiation analysis consists 
of: a), dynamic bearing force computation using the finite element and analytical gearbox 
models; b). gearbox surface velocity estimate with excitations from dynamic bearing forces 
obtained from a); and c). acoustic analysis of the gearbox with the boundary condition of 
the surface velocity calculated from b). 

The full fidelity finite element/contact mechanics model of the NASA gearbox has been 
validated. Gear transmission error is compared among commercial software and matches 
the design objective. The rolling element and gear contact properties have been validated 
against experiments in the literature. Modal analysis has been performed on the entire 
gearbox model and results agree with measurements. The simplified lumped-parameter 
model of the gearbox is suitable for parametric studies on the dynamic response. The 
analytical model has been validated against the finite element analysis. 

Cross-coupling stiffnesses are significant as they affect the gearbox dynamics. A fully- 
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Mesh Frequency, kHz 

Figure 42: Sound power of radiated gearbox noise excited by certain mesh frequency harmonics 

of bearing forces with rolling element (— ) and fluid film wave bearings ( ) at 79.09 Nm input 

torque. The excitation from 1 st to 6 th mesh frequency harmonics of bearing forces are considered 
during the computation. A weighing filter (ISO standard) is used to adjust sound pressure levels. 

populated stiffness matrix is needed to fully represent rolling element bearings. Bearing 
stiffnesses fluctuate periodically as rolling elements enter and leave the region of loaded 
rollers. These stiffness variations could excite gearbox vibration. 

This acoustic model has been validated by comparing against theoretical solutions and 
measurements taken at the NASA GRC gearbox. The transfer function method determines 
sound pressure and power for arbitrary bearings, which is exceptionally suitable for the 
gearbox optimization study to reduce noise. Noise radiation property of fluid film wave 
bearings is speed-dependent, and wave bearings are not superior to rolling element bearing 
in reducing gearbox noise. An one-way vibro-acoustic pass of the gearbox is present at the 
mesh deflection mode. Gear dynamics significantly affects the housing vibration and noise 
radiation. Housing and bearings, however, have limited effects on the gear vibration. 
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